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The  advance  of  the  front  of  hydraulic  energy  construction  into 
regions  of  Siberia  and  the  Far  East  has  generated  a need  for 
development  of  new  types  of  hydraulic  equipment,  specially 
adapted  for  use  under  the  harsh  climatic  conditions  existing 
in  these  regions.  This  approach  to  the  problem  of  the  deve- 
lopment of  hydraulic  energy  in  the  eastern  regions  of  the 
country  makes  it  possible  to  take  into  account  more  fully 
those  positive  phenomena  which  are  associated  specifically 
with  the  harshness  of  the  climate. 

The  accumulated  experience  in  planning  and  building  in 
northern  regions  shows  that  massive  dams,  erected  in  regions 
with  a harsh  climate,  develop  zones  with  a constant  negative 
temperature.  Taking  into  account  the  high  strength  charac- 
teristics and  the  favorable  deformational  characteristics  of 
concrete  in  these  zones,  we  not  only  can  increase  the  relia- 
bility of  the  structures  but  can  even  achieve  considerable 
economies. 

The  task  of  building  frozen  structures  involves  a number 
of  problems,  one  of  which  is  knowing  the  thermal  regime  of 
the  part  of  the  reservoir  which  is  in  contact  with  the  daun, 
since  it  constitutes  a warm  interface  with  the  frozen  dam 
and  foundation. 

At  the  present  time,  on  the  basis  of  studies  performed  at 
VNIIG  imeni  B.  Ye.  Vedeneyev,  Cidroproyekt  imeni  S.  Ya.  Zhuk 
emd  so  on,  so-called  "Guidelines  or  Thermal  Calculation  of 
Reservoirs"  have  been  worked  out  [1].  However,  these  do  not 
consider  the  problem  of  the  thermal  interaction  between  the 
water  and  the  structure. 
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The  goal  of  the  present  paper  is  to  provide  a general  scheme 
and  practical  method  for  calculating  the  thermal  regime  of  a ' 

reservoir,  taking  into  account  the  thermal  interaction  with 
the  dam  and  the  foundation  and  examining  the  possibility  of 
regulating  the  thermal  regime  of  the  structure. 

The  joint  consideration  of  thermal  processes  in  the  water 
and  in  the  body  of  the  structure  is  a very  complicated  problem 
in  the  theory  of  heat  exchange.  It  is  necessary  to  keep  in 
mind  that  in  the  final  amalysis,  in  order  to  determine  the 
state  of  the  structure  (including  the  frozen  parts),  it  is 
very  important  to  know  the  vertical  distribution  of  the  tem- 
perature of  the  water  near  the  structure  at  various  points  in 
time.  The  problem  of  calculating  the  temperature  of  the 
water  adjacent  to  the  structure  can  be  solved  approximately 
in  two  steps. 

In  the  first  step,  the  calculation  is  applied  to  the  tempera- 
ture of  the  water  in  the  part  of  the  reservoir  next  ti  the 
dam,  without  taking  into  account  the  thermal  interaction  with 
the  structure:  t = f(z,T). 

In  the  second  step,  the  decrease  in  the  temperature  of  the 

water  Lt  - <^(z,t),  arising  as  the  result  of  the  cooling  in-  , 

fluence  of  the  frozen  structure,  is  determined.  For  this 

purpose  it  is  necessary  to  use  the  temperature  distribution 

curves  obtained  in  the  first  step  of  the  calculation,  i 

t * f(z,T),  to  find  the  approximate  heat  exchange  intensity 

with  the  structure  Sappj-ox*  necessary,  these  results  may 

be  made  more  precise  by  performing  a second  calculation. 

First  of  all,  let  us  consider  several  aspects  of  the  thermal 
regime  of  reservoirs  which  are  important  from  the  standpoint 
of  the  thermal  aspects  of  frozen  structures:  we  will  then 
present  our  basic  views  on  the  method  of  calculating  the 
water  temperature  and  the  heat  exchange  coefficients  with 
the  atmosphere  and  the  body  of  a structure;  finally,  we  will 
deal  with  the  problem  of  regulating  the  thermal  regime  of 
structures  in  stammer. 

Analysis  of  data  from  field  observations  and  calculations 
shows  that  the  annual  thermal  cycle  of  reservoirs  contains 
five  periods  [1,2]: 

I.  spring  heating  (up  to  4”C) 

II.  summer  heating  (4'*C  and  up) 

III.  autumn  cooling  (down  to  4°C)  , 

IV.  cooling  before  ice  forms  (below  4*’C) 

V.  winter  (beneath  ice  cover) 
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Obviously,  as  far  as  the  problem  of  calculating  and  regulating 
the  regime  of  frozen  structures  is  concerned,  we  are  most 
interested  in  Periods  II  and  III  in  this  annual  thermal 
cycle,  during  which  the  water  temperature  is  highest.  Since 
the  changes  in  water  temperature  decrease  with  depth,  it  is 
advantageous  to  define  three  types  of  reservoirs  as  a 
function  of  their  depth;  shallow,  deep,  and  very  deep. 

Shallow  reservoirs  are  characterized  by  practically  homogene- 
ous temperature  throughout  their  depth.  If  no  ice  forms,  the 
change  in  the  water  temperature  as  a function  of  time  at  all 
depths  follows  the  changes  in  air  temperature,  with  a certain 
time  lag;  however,  when  the  ice  forms  the  water  temperature 
does  not  exceed  0.5  to  1°C  as  a rule. 

Deep  reservoirs  differ  from  shallow  ones  in  having  a 
temperature  which  varies  as  a function  of  depth.  During  per- 
iods of  heating,  the  water  may  be  divided  vertically  into 
three  layers;  an  upper  layer,  with  a relatively  high  tempe- 
rature; a lower  layer  with  a low  temperature;  and  an  inter- 
mediate layer,  the  "discontinuity  layer."  The  latter  is 
characterized  by  considerable  vertical  temperature  gradients. 
The  temperature  changes  as  a function  of  time  occur  in  all 
the  layers,  including  the  lower  layer.  When  the  ice  is  in 
place,  the  ten^jerature  drop  as  a f\inction  of  depth  is  2 - 3®C. 

A distinguishing  feature  of  very  deep  reservoirs  is  the  small 
cheuige  that  takes  place  in  the  water  temperature  in  the  layer 
near  the  bottom.  In  layers  located  higher  up,  the  temperature 
regime  resembles  the  regime  in  deep  reservoirs.  Using  the 
exeunple  in  Figure  1,  we  can  show  the  vertical  distribution 
and  annual  temperature  pattern  of  the  water  in  the  near-bot- 
tom part  of  the  Bratsk  Reservoir  for  1964.  As  we  see,  in 
the  layers  near  the  bottom  (beginning  at  a depth  of  30  - 40  m) 
the  water  temperature  changes  by  no  more  than  2®C  during  the 
year,  while  the  annual  variations  at  the  surface  reach  22“C. 

It  is  inqjortant  to  point  out  that  the  type  of  reservoir 
(depending  on  the  depth)  is  a function  not  only  of  the 
morphometric  depth  h,  but  also  of  certain  other  factors 
including  in  particular  the  turbulent  and  free-convective 
mixing  of  the  water,  characterized  by  the  value  of  the  coef- 
ficient of  thermal  conductivity  X. 

The  method  of  determining  the  type  of  reservoir  as  a function 
of  depth  is  presented  in  [1],  Approximately,  if  there  are  no 
data  on  field  observations  (for  example,  when  the  reservoir 
is  still  only  in  the  planning  stage)  the  type  of  reservoir  may 
be  determined  as  a function  of  depth  as  follows: 
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FIGURE  1.  Water  Tenqperature  in  the  Portion  of  the  Bratsk  Re 
servoir  near  the  dam  for  1964:  a - vertical  distribution 
(1,  April  24;  2,  May  25;  3,  July  14;  4,  August  14; 

5,  August  27;  6,  October  1;  1,  December  24);  b-  annual 
pattern  (8,  at  the  surface;  9,  at  a depth  of  10m; 

10,  at  a depth  of  20  m;  11,  at  a depth  of  45  m;  12, 
near  the  bottom) . 


’We  can  consider  that  the  reservoir  is  "shallow"  if  the  Biot 
criterion 

h 

Bi  i X/a  < 0.2, 

or  what  amounts  to  the  same  thing,  the  depth  of  the  water 

X 

h < 0 . 2 a . 

If  Bi  > 0.2,  in  other  words 

X 

h ^ 0 . 2 a 

the  reservoir  will  be  either  "very  deep"  or  "deep."  In  this 
case,  we  can  also  consider  the  value  of  the  Fourier  criterion: 

Xt 

Fo  s cph2 

When  Fo  £ 0.1  - 0.2,  corresponding  to 

r 

the  reservoir  will  be  "very  deep."  In  the  opposite  case,  the 
temperature  variations  will  reach  the  bottom  cmd  the  reservoir 
must  be  considered  simply  "deep." 

Obviously,  due  to  the  problem  of  keeping  a structure  frozen 
year  round,  the  most  vinfavorable  cases  are  those  in  which 
the  reservoir  is  of  the  shallow  type,  thermally  speaking. 

If  the  reservoir  is  a deep  one,  especially  a very  deep  one, 
the  operating  conditions  of  the  frozen  structures  will  be 
more  favorcUsle,  since  most  of  the  heat  coming  from  the  atmos- 
phere will  accumulate  in  the  upper  water  layers,  and  the  lower 
layers  will  remain  relatively  cold.  From  this  we  see  that 
the  thermal  regime  of  the  structure  can  be  regulated  by 
raising  the  colder  water  from  the  deeper  layers  to  the  sur- 
face. 

Now  let  us  move  on  to  a discussion  of  the  method  of  hydro- 
thermal  calculation. 

As  we  mentioned  earlier,  we  shall  calculate  the  water  tem- 
perature during  the  warm  period  of  the  year  (periods  II  and 
III  in  the  annual  thermal  cycle),  with  the  conditions  for  the 
hottest  summer  being  used  in  the  calculation. 
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To  calculate  the  water  temperature  in  the  part  of  a reservoir 
near  the  bottom  without  taking  into  account  the  thermal 
interactions  with  the  structure  (stage  I of  the  calculation) 
we  use  the  method  presented  in  [1],  The  original  data  are 
the  water  temperature  at  an  initial  moment  in  time  tg / compris 
ing  the  heat  exchange  with  the  atmosphere  (radiation  balance 
Sr,  heat  loss  by  evaporation  S^,  convective  heat  exchange 
Sq) , water  depth  h and  coefficient  of  thermal  conductivity  of 
the  water  (X  and  a).  It  is  convenient  to  select  tg  = 4°C  as 
the  initial  temperature  which  occurs  at  the  end  of  the  first 
pericd  of  the  cuinual  thermal  cycle;  the  corresponding  moment 
in  time  occurs  at  the  beginning  of  calculation.  As  a result, 
we  find  the  curve  for  the  change  in  water  temperature  as  a 
function  of  depth  at  any  moment  in  time  used  for  the  calcu- 
lation: t = f (z,t) . 

In  the  ser  stage  of  the  calculation,  we  must  find  the  de- 
crease i erature  At,  caused  by  the  cooling  effect  of 

the  frc 

Let  us  .jwix  up  the  thickness  of  the  water  vertically  into 
I layers  with  a constant  water  temperature  in  each  layer;  in 
other  words,  let  us  use  the  curve  which  we  have  found,  t = 

= f(z,  t)  in  the  form  of  a step  graph,  composed  of  I steps. 

To  determine  Atj^  we  first  find  the  heat  flux  values  Spp^i 
given  off  by  the  water  to  the  frozen  structure.  The  value 
of  Spei  is  determined  on  the  basis  of  cui  approximate  calcu- 
lation of  the  thermal  regime  of  the  structure,  which  here- 
after will  be  determined  more  specifically,  taking  into  ac- 
count the  water  distribution  temperature  found  in  the  first 
stage  of  the  calculations.  Two  methods  of  calculating  S 
are  possible.  ^ 


If  we  know  the  temperature  distribution  near  the  head  limit 
of  the  structure 


(1) 


where  a is  the  coefficient  of  heat  exchange;  t^  and  t^yj^  ^ 

are  the  temperatures  of  the  water  and  the  structure  in  each  of 
the  i layers,  respectively. 


If  we  know  temperatures  ti  . and  t,  . of  the  structure  at 
distances  and  X2  from  the  head  layer,  then 

• A 

‘ < 


S^m 


(2) 


where  X is  the  coefficient  of  thermal  conductivity  of  the 
structure. 
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-hen,  knowing  the  heat  flux  Sp^  £ » we  determine  the  decrease 

in  water  temperature  At.  at  various  distances  from  the  struc- 
ture at  individual  momeAts  in  time.  This  problem  corresponds 
to  the  conditions  in  problem  No.  15  in  [1].  The  temperature 
drop  in  each  layer  will  be 

(3) 

where  x is  the  distance  from  the  dam.  To  determine  the  tem- 
perature parameter  0 we  use  a calculation  graph  (Figure  2) 
where  the  original  argximent  is  Fo  = ax  . 


Basically,  for  each  of  the  layers,  we  find  the  temperature 
of  the  water  ty,  i = - At-  at  individual  moments  in  time 

at  various  distlinces  from  the  dam,  carrying  out  our  calcula- 
tions to  a distance  x , at  which  At.  = 0,  i.e., 

approx,!  i 
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FIGURE  2.  Calculation  Graph 
to  Determine  Water  Temperature 


The  method  of  calculating  the  influence  of  the  frozen 
foundation  on  the  water  temperature  is  theoretically  the  same 
as  that  set  forth  above.  The  decrease  in  water  temperature 
produced  by  heat  exchange  with  the  frozen  foundation  will  be 


(4) 


The  temperature  parameter  9 is  determined  by  the  graph 
(Figure  2)  , where  Fo  5 at  . The  values  of  S.  are  indi- 

(h-z)  2 

cated  approximately,  taking  into  account  the  hydrothermal 
calculation  of  the  water  temperature  near  the  bottom  computed 
in  the  first  step. 

The  method  set  forth  for  calculating  the  water  temperature 
does  not  take  into  account  the  water  temperature  variations 
from  day  to  day.  However,  in  regions  with  a continental 
climate  during  the  warm  period  of  the  year,  day-to-day 
temperature  variations  are  consideraQale  and  may  turn  out  to 
ha/e  a considerable  influence  upon  the  thermal  regime  of 
frozen  structures. 

To  determine  the  aunplitude  of  the  diurnal  variations  in 
water  temperature,  this  diurnal  pattern  f.or  air  temperature 
must  be  represented  as  em  equivalent  cosine  curve  with  a 

period  tn  * 24  hours;  0 = 9 + T„  cos  2tt— ^ . 

^ Cp  0 To 

Then  [3]; 
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Where  T is  the  amplitude  of  the  water  temperature  variations 
at  depth  z;  To  is  the  amplitude  of  the  temperature  varia- 
tions of  the  air;  0 is  the  running  air  temperature:  ©cp 
is  the  average  daily  air  temperature;  Bi^FOg  e ; 

Bi  = «a2  is  the  Biot  criterion;  Fo-  = atp  is  the  Fourier 
X z2 

criterion;  oia  is  the  coefficient  of  heat  exchange  at  the 
water-air  interface;  z is  the  distance  from  the  surface  of 
the  water. 


FIGURE  3.  Calculation  Graph  to  Determine  Daily 
Variations  in  Water  Temperature. 


Using  this  formula,  a calculation  graph  was  plotted 
(Figure  3)  which  can  be  used  to  find  the  amplitude  of  the 
water  temperature  variations  at  a depth  z (at  known  values 
of  Bi^FoQ  and  Foq): 


I* 


(6) 


or  the  thickness  of  the  layer  at  whose  lower  limit  the  ampli- 
tude of  the  water- temperature  variations  constitutes  a defi- 
nite part  of  the  amplitude  of  the  air-temperature  variations 
(at  known  values  of  Bi2FoA  and  A = T/Tq) s 

(7) 

As  the  air  temperature,  we  must  use  the  "effective  tempera- 
ture : " 


r S, 


(8) 


Now  let  us  discuss  the  question  of  the  labelling  of  the 
heat-exchange  coefficients  at  the  water-air  boundaries 
(ug)  and  the  water-structure  boundaries  (a) . 

To  calculate  the  heat-exchange  coefficient  between  the  water 
and  the  atmosphere,  we  can  recommend  the  formula  of 
A.P.  Braslavskiy  et  al.,  presented  in  [1]: 

- 2iC(l  ; 0,8Br  r f {i4)\,  KKOJlljfl-l.  (9) 

Here  W is  the  wind  speed,  m/sec;  f(/i0)  is  a function  which 
takes  into  account  the  influence  of  the  stratification  of  the 
air  near  the  water  layer  (the  values  of  f(A9)  are  given  in 
[1]). 

The  question  of  the  significance  of  the  heat  exchange  co- 
efficient between  the  water  amd  the  structure  is  one  that 
has  been  little  investigated.  In  the  general  case,  heat 
exchange  may  take  place  due  to  both  forced  and  free  convec- 
tion. To  determine  the  heat  loss  coefficient  with  forced 
convection  it  is  necessary  to  know  the  rate  at  which  water 
flows  around  the  structure.  An  analysis  of  materials  on  the 
hydraulic  regime  of  existing  reservoirs  shows  that  in  the 
majority  of  reservoirs  the  flowrates  are  very  insignificant 
in  the  part  near  the  dam.  An  increase  in  velocity  occurs  only 
in  a zone  located  near  the  water-collecting  structures 
(spillways,  bottom  openings,  penstocks,  etc.).  The  size  of 
this  zone  is  small  as  a rule  by  con^arison  with  the  total 
area  of  the  head  front  structure.  This  indicates  that  the 
principal  role  in  heat  exchange  between  the  water  and  the 
structure  is  played  by  free  convection. 
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The  value  of  the  coefficient  of  heat  loss  a may  be  determined 
using  the  formulas  of  M.A.  Mikheyev,  which  relate  to  the  case 
of  free-convective  heat  exchange  around  a wall  [41: 

''»ith  Ra  = 1 X 10-3  to  5 X 102  Nu  = 1.18  Ra^/S 

With  Ra  = 5 X 102  to  2 X lO"^  Nu  = 0.54  Ra^/^  (10) 

With  Ra  = 2 X 10*^  to  1 X 10^^  Nu  = 0.135  Ral/3 

3 

Here  Nu  i ^ is  the  Nusselt  criterion;  Ra  = gBAtd  ir. 

the  Rayleigh  criterion;  g is  the  acceleration  due  to  gravity; 

3 is  the  coefficient  of  thermal  expansion  of  water; 

At  = ti  - t d is  a critical  parauneter  equal  to  the 
thiclcness  or  the  I-th  layer  of  water;  V(j),X4i,  and  a^  are 
physical  coefficients  of  kinematic  viscosity,  heat^and 
temperature  elasticity  of  water. 

We  mentioned  earlier  that  one  question  of  great  practical 
significamce  is  the  regulation  of  the  thermal  regime  of  fro- 
zen structures  during  the  summer  period,  owing  to  the 
fact  that  colder  water  rises  to  the  surface  from  the  deeper 
and  near-bottom  layers.  The  possibility  of  hydrothermal 
effects  upon  the  thermal  state  of  the  structure  arises  if 
the  reservoir  is  of  the  deep  or  very  deep  type.  In  such 
reservoirs  the  temperature  near  the  bottom  may  be  15  - 20 ’C 
cooler  them  the  surface  layers  (Figure  1) . 

To  lift  the  bottom  water  toward  the  surface,  pumps  may  be 
used,  as  well  as  flow  converters  and  pneumatic  units.  The 
effectiveness  of  these  measures  can  be  increased  by  installing 
vertically,  along  the  front  of  the  structure  at  a distance 
of  0.5  to  1 meter,  u lightweight  structure,  or  "curtain" 

(for  exanple,  a polyethylene  film  mounted  on  the  body) , which 
does  not  reach  the  bottom  of  the  reservoir.  The  upper  edge 
of  this  curtain  must  be  slightly  depressed  below  the  surface 
of  the  water.  This  curtain  considerably  decreases  the  water 
and  heat  exchange  between  the  cavity  it  forms  emd  the  rest  of 
the  reservoir.  An  even  more  effective  way  is  to  use  a 
double  curtain,  which  decreases  heat  loss  through  both  con- 
vection emd  thermal  conductivity. 

In  conclusion,  we  should  like  to  point  out  that  preliminary 
calculations  using  the  method  outlined  above  have  shown  that 
in  the  constant  temperature  area  the  intensity  of  the  heat 
exchange  between  the  reservoir  amd  the  structure  may  amount 
to  0.5  kW/m2.  In  the  variable  teirperature  zone  (down  to 
10  - 20  m)  the  intensity  of  the  heat  exchange  is  usually 
2 to  3 orders  greater.  Therefore,  it  must  be  considered 
necessary,  when  it  is  desiradsle  to  have  water  temperatures 
which  constitute  the  boundary  condition  in  calculating  the  de- 
sign of  structures,  to  take  into  account  the  heat  exchange 

at  the  water-structure  interface. 
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